Modelling gravity on a hyper-cubic lattice 
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We present an elegant and simple dynamical model of symmetric, non-degenerate n x n matrices 
of fixed signature denned on a n-dimensional hyper-cubic lattice with nearest-neighbor interactions. 
We show how this model is related to General Relativity, and discuss multiple ways in which it can 
be useful for studying gravity, both classical and quantum. In particular, we show that the dynamics 
of the model when all matrices are close to the identity corresponds exactly to a finite-difference 
discretization of weak-field gravity in harmonic gauge. We also show that the action which defines 
the full dynamics of the model corresponds to the Einstein-Hilbert action to leading order in the 
lattice spacing, and use this observation to define a lattice analogue of the Ricci scalar and Einstein 
tensor. Finally, we perform a mean-field analysis of the statistical mechanics of this model. 
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Introduction: Lattice models of gravity are typically 
defined using some discretization which is simultaneously 
coordinate and background independent. The most pop- 
ular such discretization is the Regge calculus [H, [3 .used 
to study both classical general relativity (GR) [3j, |4| , as 
well as to build models of quantum gravity based on dy- 
namical simplices (B- 17|. Another discretization, based 
on a topological lattice, is [ijj], (see also [lj|), however 
this discretization is designed to preserve some form of 
diffeomorphism invariance, and so is also both coordinate 
and background independent. 

In this Letter we shall instead present a discrete model 
for gravity which is defined on a regular, in fact hyper- 
cubic, coordinate lattice. The implied background struc- 
ture may be anathema to GR purists, however we will 
argue that this is still a useful thing to do, and can be 
usefully utilized to study GR. 

That having a preferred background might not be en- 
tirely implausible may be inferred from a number of 
observations: (1) Many interesting spacetimes can be 
put into Kerr-Schild form, which has a natural back- 
ground [20j . (2) Many interesting spacetimes can be put 
into Painleve-Gullstrand 2l| and/or de Donder (non- 
linear harmonic) form (22[, both of which possess natu- 
ral background metrics. (3) Many interesting spacetimes 
can be put into "relativistic acoustic" form, based on the 
"analogue spacetime" programme, for which a natural 

(4) Physically in- 
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background metric again exists 
teresting black holes can be put into horizon-penetrating 
coordinates, for which the metric components are finite at 
the horizon; the presence of a horizon does not necessar- 
ily imply "infinite deviations" from some assumed back- 
ground metric [21, 24 1. (5) More exotically, recent spec- 
ulations on ghost-free massive gravitons are most nat- 
urally phrased in terms of a combination of foreground 
and background metric 25]. In view of the above, we 
are willing to at least entertain the notion of background 
structure, to see how far we can get. 



Lattice action: Consider a n-dimensional hyper-cubic 
lattice which has defined, at each site i, an n x n sym- 
metric, non-degenerate matrix l g, which is physically to 
be interpreted as the metric. Unless otherwise stated, 
we will assume in this Letter that the matrix is positive 
definite, (and hence is a model for a Euclidean-signature 
Riemannian geometry), however the model can easily be 
generalized to matrices of any fixed signature. The dy- 
namics of the model is described by a particularly simple 
and elegant action defined as a sum over nearest neigh- 
bor pairs. Let us first define an "average" metric linking 
the sites i and j 
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and then set 



s = -i n ( s /L P y 
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Here (ij) denotes the link joining nearest neighbor 
sites i and j. The lattice spacing is s, and the n- 
dimensional Planck length is Lp, while is some conve- 
nient dimension-dependent normalizing constant. Note 
that this action has a symmetry under both rigid (global) 
SO(n) transformations % g — > O t % g O, and under parity. 
Weak field: Consider small fluctuations 

l g = l + \ (3) 

To quadratic order in h we have 

V^W) = 1 + \ tr [ l h + j h] + ^ tr [Vi + %] 2 

-j-tr[Q h+ j h f]+0(h 3 ), (4) 
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and 
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- v/dctM = -+- tr [%]+_ tr [%] --tr [ l h 2 ]+0(h 3 ). 



Thus to quadratic order in h we have 
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ftr [ (% - l hf] - itr %] 2 ^) + 0(/i 3 ). 

(6) 

One might worry that the matrices were only expanded to 
linear order; however repeating the calculation and keep- 
ing the quadratic terms gives the same result. Taking 
the lattice spacing to zero the finite differences become 
derivatives, and the action (up to an arbitrary multiplica- 
tive constant) is given by 

s cx J d n x (d a h^ srhr - i d»h% d^J) +o(h 3 ). 

(7) 

This is precisely the action for linearized GR in (lin- 
earized) harmonic gauge 



d v Ku - I d,h\ = 0. 



(8) 



Thus we see that the lattice action (|2J) can be used to 
model weak field GR provided the very commonly used 
(linearized) harmonic coordinates are adopted. 

Continuum limit: Naively, it might seem just a lucky 
coincidence that the coefficients in the linearized expan- 
sion are precisely those needed to recover the GR terms 
of weak-field gravity. However, we will show that the 
discrete lattice action ([2]) has a continuum limit, which 
is given to leading order by the Einstein-Hilbert action. 
Let us use permutation symmetry to first rewrite the ac- 
tion as 



S cx £ { VfctWv) - TdetM} 



We can explicitly factor out a yj det( l <?) to obtain 



(9) 
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Here we define the site-specific contribution to the action 
for site i in terms of a sum over its nearest neighbours 

ft (x £ (2-«/ 2 v^fl+ira - l) • (11) 

Let us now work on a continuum manifold with metric 
9tiu{x) and choose a Riemann normal coordinate system 
— such that site i is taken to be the origin, and the metric 
at site i is l g^ v — 6^. Let £^ denote the unit vector 
pointing from site i to site j. The coordinate system 
is such that the geodesies generated by these coordinate 
unit vectors are straight lines: 



v x tl (\) = i?j A. 



(12) 



Then in the immediate neighborhood of the origin we 
can construct the vertices of the hyper-cubic lattice such 
that the nearest neighbors are connected by geodesies of 



(5) length s and have coordinate locations J a; M 



M <? A 

standard result for Riemann normal coordinates is that 
to quadratic order 



V = S/u^aO = <W - o ?WjS 1% 1% s 2 + 0(s 3 ), (13) 



But then 

{[W 1 = *»* - I K»a»0 t% 4 S 2 + 0(S 3 ), (14) 

so 

2-"/Vdct(I+fe]- 1 [Jg]) = 1-1 n aP 1% 4 s 2 +0{s 3 ). 

(15) 

Hence the site-specific contribution to the action is 

*S cx K a p 1% 4 s 2 + 0(s 3 ). (16) 

Using the easily deduced result 

E^4 = 2 ^ ( 17 ) 

we see that to quadratic order in the lattice spacing 
(which is also the geodesic distance) the continuum ana- 
logue of equation (|lip is given by 



l S otKs 2 +<D(s 3 ). 



(18) 



Thus we see that our lattice action corresponds to the 
Einstein-Hilbert action, to leading order in the lattice 
spacing. 

Lattice Ricci scalar: Choosing a suitable normaliza- 
tion we define 

lR = - 6 E ( 2 ~ " /2 \/ditOf+feP M) l) • (19) 

Then trivially adapting the discussion above 

l R = Ks 2 +0(s 3 ). (20) 

We see that to lowest nontrivial order in the lattice spac- 
ing, the discrete quantity l R exactly matches its contin- 
uum analogue 1Z. 

Strong field EOM: The discrete version of the Ein- 
stein tensor is easily obtained by computing 
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SS 
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(21) 



We find 



|4£ (v^W) C'v)- 1 

Lyj r-(ij) \ 
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(22) 
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After picking a suitable normalization, we set 

rv= - o > , — — '- r-r 1 -^)- 1 




(23) 



But now (again adopting Riemann normal coordinates at 
the site i, so l g — > I) we have already seen (to quadratic 
order) the equivalent of 

1 

6 

and 



[ ij 9] ll = <W + £ Kvcp 1% 1% s 2 + 0(s 3 ), 



det [Vg]=l--K a0 t%l%s 2 + Q{s 3 ) 
whence, to quadratic order 



(24) 
(25) 



^ - 2 



+ 0{s 3 ). 
(26) 

So, summing over the nearest neighbour sites j, the dis- 
crete Einstein tensor is related to the continuum Einstein 
tensor by 



IS1 

"(11/ 



s 2 +0(s 3 )=G^ s 2 +0(s 3 



(27) 

External stress- energy tensor: As is usual in lattice 
models, we can also add an external current *J to probe 
the dynamics. The most natural object to add is 



(28) 



The external current l J is interpretable in terms of the 
discrete stress-energy tensor l T via 

'T = - 1 = V + I fg}- 1 trPg % (29) 

x/detM %] 2 LVJ Ly 

The strong field discrete Einstein equations in the pres- 
ence of external stress-energy are then quite simply 



(30) 



So formally at least, the discrete lattice model contains 
all the correct ingredients for adequately dealing with 
large swathes of standard GR. This procedure clearly 
generalizes to placing some matter model on the lattice. 

Mean-field analysis: In addition to using this lattice 
model to study classical GR, it can also be used as a 
discrete model for studying quantum gravity. A first step 
in this direction is to perform a mean-field analysis of 
the action ^ which we now rewrite in terms of the site- 
specific form as 



S=-a n J2 y/fetft) (31) 

i 

(2-/ 2 Vdet(fTfeP M) l) » 



r-{ij) 



where a r , 



Cn(s/Lj 



\n-2 



is the constant appearing in 



equation ([2]). This action is translationally invariant 
and thus we take the mean-field ansatz — assuming the 
physics is dominated by some translation invariant aver- 
age M — (g), plus small fluctuations 



i g = M + 8[ i g). 



(32) 



This allows us to replace 3 g with M in the coupling term 
of the action. We find 



fJ\ 



1 



-> 2n [2- n / 2 v/det (1+ fe]" 1 M) - 1 
The total action then becomes 

Smf = -«« 2n 2_j V det (1?) 

i 

x (V "/ 2 Vdet (1+ fe]" 1 M) - l) 
Thus the mean field partition function 



(33) 
(34) 
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is given by 

Z [ / dffe 7^(2-/V«ict(I +ff -M)-l) 



TV 



(35) 



(36) 



where 7 = 2na n (3. Now performing a change variables 
g = \fM g\J~M \ which has Jacobian det J = (detM)™, 
after dropping the tilde the integral which determines 
the mean field partition function is given by 



/ 



dg e 



(37) 



Thus we see that our mean field analysis results in a ran- 
dom matrix model which is invariant under 0{n) trans- 
formations. Because the matrix model has this symmetry 
group we can perform the diagonalization g — T A O 
where A is the matrix of eigenvalues, all of which are by 
assumption positive. The Jacobian of this transforma- 
tion is given by [H, ^ 



dg = dO\{d\ A U \\a-\b\ 



(38) 



A<B 



Thus, after integrating over the orthogonal group, the in- 
tegral appearing in the mean field theory partition func- 
tion becomes 



[l[dX A J] \\ A -\ B \e-tV^WvM, 

J A A ^ D 



A A<1 

where the function V(A) is given by 



1 + Xa 



(39) 



(40) 
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After a bit of work we obtain the (useful but sub-optimal) 
constraints 

n\ / Ai-n^ <F w < n^ w 

A A <1 A a >1 A 

Only for the trivial case n = 1 (a one-dimensional chain) 
is the function V(A) bounded from below; for n > 2 there 
are directions in eigenvalue space where V(A) becomes 
arbitrarily negative. Thus the mean field truncation of 
the lattice action gives a random matrix model exhibiting 
pathology similar to that of the Einstein-Hilbert action. 
Further study of this model, using the techniques of ran- 
dom matrix theory, may shed light on how to deal with 
this feature. 

Discussion: We have seen that with the particularly 
simple and elegant discrete action ((2]) one can successfully 
encode a very large fraction of standard GR. The action 
is gauge fixed, with only rigid (global) rotations and par- 
ity inversions as symmetries, and seems automatically to 
be in the de Donder (nonlinear harmonic) gauge; this is a 
feature, rather than a problem — the hyper-cubic lattice 
and the de Donder gauge seem related at some deep level. 
Presumably there is some more general gauge-invariant 
action of which this is gauge-fixed version. We note also 
that this model does not exhibit the active diffeomor- 
phism symmetry of the continuum theory and thus there 
are n(n + l)/2 dynamical degrees of freedom per lattice 
site; this is standard for discretisations of gravity [jjj ]. 
Of course gauge fixing is not a problem per se fjjj , since 
to make physical predictions one ultimately has to do so 
anyway. What is perhaps a little surprising is just how 
far one can get with such a simple and elegant discrete 
action. 
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